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(54) A method for the supression of multiple reflections from marine selmsfc data 



(57) A method for the suppression of multiple reflec- 
tions using a Kirchhoff algorithm ia daaoribod- The 
method allows the simulation of the free-surface multi- 
ples out of a zero-offset or input section. The method 
steps comprise: i} simulation of multiple reflections by 
means of a Kircnnoff-type summation applied to a 



pseudo zero-offset section; li) application of an adaptive 
filter to adjust the simulated multiple and Hi) subtraction 
of the so-obtained adjusted simulated multiple from the 
input data. 



FIGURE 2 
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Description 

FIELD OF THE INVENTION 

5 [0001J The present invention relates to a method for the suppression of multiple reflections from marine seismic 
data, more specifically, to a method for the near offset range suppression of multiple reflections where the simulation of 
multiple reflections by means of a Kirchhoftoype algorithm is applied to a common near-offset section. As a result, sig- 
nificant Improvement in the seismic Image Is achieved. 

10 BACKGROUND INFORMATION 

[0002] Routine interpretation of processed seismic data is mainly focused on understanding primary reflections. 
Multiples, however, also occur in ail seismic reflection data, acquired both in land and marine environments. Although 
multiples may bring useful data on the subsurface medium, present processing technology is still not able to make good 
is use of this information* On the contrary, multiple arrivals are seen as undesirable events. 

[0003] If not properly attenuated or removed, multiples can interfere with primary reflections, leading to difficulties 
in the estimation of velocities and amplitudes. In many land data cases, multiple attenuation can be adequately 
achieved by deconvoiution and CMP stacking. 

[0004] Because of the presence of strong reflections on the water-atr contact (free surface), in particular sea-bot- 

20 torn multiple reflections, marine data most often require additional processing. 

. [0005] in general, methods of multiple attenuation or suppression can be divided into two categories, namely (a) 
filtering - based on some characteristics of the multiple energy involved and (b) prediction - simulation of multiples and 
subtraction from the original data. The welJ known Radon and f-k transform algorithm*, generally applied in the CMP 
domain, belong to the first group. An internal mute on near offsets has to be applied, as a discrimination between prh 

25 maries and multiples is more difficult in this range. Algorithms of the second category are baaed on wave-propagation 
approaches such as forward modeling or inverse scattering. These are generally expensive and require, moreover, a 
good sstimate of the source signature for useful results. 

[0006] The increasing use of seismic data in reservoir characterization, together with the large number of explora- 
tory prospects on targets of stratigraphic nature, poses as a very desirable aim any improvement in resolution that can 
30 be achieved by seismic processing. In this context accurate elimination of multiple reflections in near-offset range 
within a target region may be crucial to correctly image, say. a fine stratigraphic sequence or to derive reliable seismic 
attributes of a reservoir. 

[0007] In the case of deep-water data, suppression of first-order sea-bottom multiple and peg-leg reflect Jons (i.e., 
which reverberate only once in the water layer) Is of great importance. These rather strong multiples may have the same 

xs traveitime as the primary reflections of target reflectors. 

[0008] US patent4.665.5t0 teaches a method of attenuating multiple reflections in a conventional seismic process- 
ing flow. Predictive Deconvoiution After Stack is used to reduce the multiple interference on stack velocity analysis. The 
Deconvoiution After Stack emphasizes the primary reflections with their stack velocity ranges, and this information can 
aid in selecting a good multiple stack velocity for filtering in F-K domain. 

40 [0009] US patent 5.729.506 teaches a computationally economical method for applying multi-dimensional multiple- 
reflection attenuation to a marine seismic-signal data set comprised of a plurality of common shot sections that have 
surface-multiple wavefields embedded therein. A common shot section is chosen from among the plurality of common 
shot sections. The predicted surface-multiple wavefield nearest to the chosen common shot section is adaptively fil- 
tered to match the predicted surface-multiple wavefield to the surface-multiple wavefield embedded in the chosen con> 

45 mon shot section. The matched predicted surface multiple wavefield is subtracted from the embedded surface multiple 
wavefield to provide a multiple-free common shot section. 

[0010] US patent 5,051 ,960 teaches a method of removing multiple reflection events from seismic records. A shal- 
low primary reflection event is used as a model for a deeper primary, and is distorted rsspon&ive to the predicted static 
correction to yield a model for the multiple. The method is applied to a common depth point (CDP) record, in which the 
50 common-midpoint contribution of the multiple to aaoh CDP record can bo assumed to be equal. Th« technique taught 
in said US patent is applied soleiy to removal of multiples which arise from reflection at the surface during land explo- 
ration. 

[001 1] Rubral. P. at al in "A unified approach to 3-D seismic reflection imaging. Part I: Basic Concepts*, Geophysics 
vol. 61 n° 3 p.742-758 (1066) state that the basic idea of the unified approach is to cascade Kirebhoff migration and 
55 demigration integrate to solve a variety of seismic imaging problems. The application of both integrals in sequence leads 
to weighted summation along certain stacking surfaces. This kind of algorithm is known in the technique as a Kirchhoff- 
type algorithm. The present invention teaches one possible application of the principles discussed in said article. 
[0012] Tygd, M. ot al. in "A unified approach to 3-D seismic reflection Imaging. Part II: Theory", Geopnyslcs vol. 61 
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n* 3 p. 759-775 (1996) discuss the theoretical principles stated in Part I of the article. 

[0013] Wegiein, A. B, in "Multiple attenuation: an overview of recent advances and the road ahead (1999)", The 
Leading Edge, January 1999, pages 40-44, presents an overview of the current state of multiple attenuation as well as 
recent advances in the field. Filtering as well as prediction methods such aa wavefteld prediction and subtraction, wave- 

5 field extrapolation, free-surface multiple elimination, feedback and inverse-scattering methods are discussed. 

[0014] Filpo Ferreira da Silva, E. and Tygel, M., In "Deep-water multiple suppression in the near-offset range**- The 
Leading Edge, January 1999, pages 81-84, present a method to suppress first-order and peg-leg multiples in deep- 
water marine data. The main steps are: t) simulation of multiples py means of a weighted Kirchhoff-type summation 
applied to a stacked section; ii) application of an adaptive filter to adjust the simulated multiple: and iii) subtraction of 

w the adjusted multiple from the original data. Advantageously, the technique has no limitations on short offsets and no 
requirement to estimate source wavelets and reflection coefficients. 

[0015] Spitz, S., in "Pattern recognition, spatial predictability, and subtraction of multiple events", The Leading 
Edge, January 1999, page 55 to 58, discusses the extraction of simulated or modeled multiples from the original data 
set, which corresponds to step iii) of the method described by Filpo Ferreira da Silva and Tygel M. In the previous par- 
ts agraph. En this article Spitz shows that the subtraction methodology is mainly concerned with the shapes of the events. 
[0016] Filpo Ferreira da Silva. E. and Tygel, M. ( in "Near-offset multiple suppression", Society of Exploration Geo- 
physicists (SEG), Sixty-ninth Annual Meeting. Tulsa. USA. October 31 -November S, 1999, describe a simulation 
method to effect a transformation of a given common near-offset or stacked input section (taken as an approximation 
of a zsro-offset section and containing all primaries and multiples), into a corresponding output section that solely con- 
20 sists of the multiples to be suppressed. 

SUMMARY OF THE INVENTION 

[0017] Broadly, the main steps of the application of the method are j) simulation of multiple reflections by means of 
za a Kirchhoff-type summation applied to a pseud o zero-offset section. As explained hereinbefore, this Is an Input common 
near-offset of stacked section that is taken as an approximation of a zero-offset section; ii) application of an adaptive 
filter to adjust the simulated multiple reflections and 111) subtraction of tne obtained adjusted simulated multiple reflection 
from the original data. 

[0018] In the present invention, steps 0) and iii) are carried out according to state-of-the-art principles. 
30 [0019] Thus, the present invention provides a method for the suppression of multiples with the aid of a Kirch hoff 
algorithm. 

[0020] The present invention provides further a method where the suppression of multiples may be applied to the 
near-offset range. 

[0021] The present invention provides stilt a method for the suppression of multiples which may be applied to mut- 
3£ tiples generated by sea bottoms of any dips. 

[0022] The present invention is based on theoretical principles which lead to a simulation of multiples in a 3-D envi- 
ronment even when using streamer data, 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 

[0023] 

FIGURE 1A to 1D illustrates the method of the present invention, in graphical form. 

FIGURE 2 illustrates the method of the present invention as a flow sheet, with each corresponding step of Figure 
45 1 being indicated. 

FIGURE 3 illustrates the Kir oh hoff method of multiple simulation. 

FIGURE 4 illustrates an embodiment of the method of the invention for the suppression of multiples, 
FIGURE 5 illustrates another embodiment of the method of the invention for the suppression of multiples. 
FIGURE 8 illustrates a further embodiment of the method of the invention for the suppression of multiples. 
so FIGURE 7 illustrates a still further embodiment of the method of the invention for the suppression of multiples. 

FIGURE 8 is an example of a seismic section where a target horizon is affected by a multiple reflection. 
FIGURE 9 is a zoom view of the target horizon of FIGURE 7. 

FIGURE 10 is FIGURE 9 after applying the present method for the suppression of multiples. 
figure n Is an illustration of the method of the Invention on a stack section. - 
55 FIGURE 12 A is a 2oom view of FIGURE 1 1 before applying the method of the invention while FIGURE 12B is a 
zoom view of FIGURE 11 after applying the method of the invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

(0024] Broadly, the method of the present invention teaches an accurate and efficient scheme to simulate free-sur- 
face multiples out of a zero-offset (or input) section (inoluding primaries and multiples), A reasonable approximation of 

s the original zero-offset section can be a short-offset or a CMP-stacked section. 

[0025] Thus, the main steps of the inventive method are i) simulation of multiple reflections by means of a Kirchhoff- 
type summation applied to a pseudo zero-offset section; ii) appiication of an adaptive filter to adjust the simulated mul- 
tiple and iii) subtraction of the obtained adjusted multiple reflection from the pseudo zero offset section. FIGURE 1A to 
1 D and FIGURE 2 illustrate said main steps. 

10 [0026] in the present specification, the expression "pseudo zero off section" means a section that contains all pri- 
maries and multiples and can serve as a good approximation of a zero offset section. For example, a common offset 
section transformed Into a zero offset section or a stack section. 

[0027] The application of the present method leads to a weighted Kirchhoff-type stack which, applied to the pseudo 
zero-offset section, performs the task of transforming an Input zero-offset section, that comprises primaries and multf- 
is pies, into an output simulated section In which ail original primary reflections have been transformed into their corre- 
sponding, free-surface multiples. 

[0028] It should be noted that the weight function that is to be applied in the present Kirchhoff-type method can be 
chosen so that the resulting amplitudes of the multiples welt approximate their observed values. The use of this weight 
function should provide improved initial amplitudes prior to the application of the filtering procedures. 

20 [0029] Thus, the simulation method described in the present application is designed to effect a transformation of a 
given common near-offset or stacked input section (the input pseudo zero-offset section), into a corresponding output 
section that solely consists of the multiples to be suppressed, As stated before, the input pseudo zero-offset section is 
taken as an approximation of a zero-offset section and containing all primaries and multiples. 
[0030] According to the illustration of FIGURE 3. in a Kirchhoff-type stacking algorithm, for example a Kirchhoff sim- 

25 ulation of multiple reflections, each sample of the output section (that is. the simulated multiple section) is the result of 
a sum of samples collected along a uniquely specified curve in the input section (the pseudo zero offset section). This 
curve is generically called the stacking line. The input samples along the stacking line are multiplied by certain weights, 
also determined by the sample in the output section, where the resulting sum is to be placed. The terminology "weight 
function" is employed to refer to the choice of multipliers that are used in the summation process, which in turn is called 

30 a weighted Kirchhoff stack. In symbols, the Kirchhoff-type transformation can be written 

Wherein the following notation has been used 

40 fo.fo) is the midpoint and traveltime coordinates of the output point in the simulated section of multiples; 

^o(Su.'o) Is the amplitude resulting from the application of the algorithm to be placed at the point ($o.fo)m the out- 
put section; 

(^,0 is the midpoint and traveWme coordinates of the points in the input pseudo zero offset section; 

is the amplitude at the point ft,f) in the input section; 
45 WfcSo.fo) & the weight function; 
t&£o.'o) is the stacking line; 

O|f[WR.0] is the Hilbert transform of the temporal half derivative of Lf(£.f). 

[0031] The crucial step in any Kirchhoff-type summation algorithm is the construction of stacking lines and weight 
so functions. This also applies to the present simulation scheme, with the exception that weiahts are not so crucial due to 
the further application of an adaptive filter to adjust amplitudes and phases. The use of adaptive techniques to match 
discrepancies between predicted and observed multiples is a common procedure in schemes of multiple attenuation. 
[0032] For a fixed point (5o. fo) in the output section, the expression of the corresponding stacking line drawn in the 
input section can be shown to be as stated in equations (2), (3) and (4) below: 



S3 



(2) 
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where 
in which 

♦ -^-'(s^b) (4) 



v being the p-wave velocity propagation. 

[0033] The above formula is an exact traveltime description of the firs border sea-bottom multiple in the zero-offset 
15 situation. It provides also a good approximation of the traveltime* of corresponding near sea-bottom peg-leg multiples. 
[0034] The simulation of the correct amplitude and pulse shape of the multiples to be suppressed can be theoreti- 
cally achieved by the Introduction of weights in the Kirchhoff-type stacking integral of equation (1). Those weights 
depend, however, on the reflection coefficients and incident angles of the multiples to be suppressed. As these quanti- 
ties ore not known, it ia much more practical to effect the necessary amplitude adjustments upon the uao of adaptive 
20 filtering techniques. This procedure is also helpful in handling peg-leg multiples. 

[0035] The simulation of a multiple reflection can be achieved by a cascaded application of Migration and Demigra- 
tion processes. First, the pseudo zero-offset section is migrated to the depth domain- Next, the obtained depth section 
is demigrated back to time domain so as to obtain the simulated multiple section. 
[0036] Thus, the zero-offset Kirchhoff migration is performed by the integral as in equation (5): 



Wherein the following notation has been used: 



(5) 



A M is the migration aperture: 

35 (x. 2) is the spatial coordinates of the output point in the depth migrated section; 

V{x. z) is the amplitude resulting from the application of the depth migration algorithm to be placed at the point 

(x, z) in the output migrated section; 

is the midpoint and traveltime coordinates of the points in the input pseudo zero offset section: 
U{%,t) is the amplitude at the point {£,/) in the input section; 

40 Wm (£; x. z) is the weight function of the migration algorithm; 
x M (£; x, z) is the stacking line of the migration algorithm; 

Df{D(£,f)] io the H Libert transform of the temporal half derivative of U(q,t). 



The stacking line x w ft; x,z) in the migration algorithm is given by equation (6) 

t w (^^(x.5)^z z (6) 



so [0037] The weight function (£; x t z) in the migration algorithm te taken as unitary In this procedure. 

[0038] The zero offset Kirchhoff demigration algorithm used to obtain the multiple section is performed by the inte- 
gral (7) 



55 
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10 



wherein 



A D is the demigration aperture; 

(£>,fo) i$ the midpoint and traveJtime coordinates of the output point in the demigrated section, which corre- 

sponds to the simulated multiple section; 

U(£o.'t>) is the amplitude resulting from the application of the demigration algorithm to be placed at the point 
(qo.to) in the output demigrated section; 

(x t z) are the spatial coordinates of the points in the input depth migrated section; 

V[x t z) is the amplitude at the point (x*z)in the input depth migrated section: 

WjD(x.t^o,/o) ' s the weight function of the demigration algorithm; 

to(*.£o*fo) is the stacking line of the demigration atoorithm: 
20 0 2 [V{x,2)) is the Hilbert transform of the depth half derivative of VT(x.z). 

[0039] The stacking line tofx^o.'o) >n the demigration algorithm is given by equation (8) 

™ t fx* t m 

[004Q] The weight function Wp(x t c^j 0 ) in the demigration algorithm is also taken as unitary in this procedure. 
{0041J Aa illustrated in Figures 4, S t 6 and 7, an ensemble of common offset sections is made up of several sec- 
30 trons. each section being made up of several seismic traces registered with the same distance between source and 
receiver. As we J I known by the experts, the distance between source and receiver Is designed by the expression "off- 
set", such expression being used hereinbelow. 

[0042] The several sections are organized In increasing order of offset or distance, from the shortest to the largest 
offset. 

35 [0043] In one embodiment, as Illustrated In the flow sheet or FIGURE 4, the method of the invention for the suppres- 
sion of multiple reflections from a marine seismic data set having registered primary and multiple reflections comprises 
the following steps: 

a) sorting said input data set so as to obtain an ensemble of common offset sections; 
40 b) transforming each of said common-offset section into a zero offset section, so as to obtain an ensemble of zero- 
offset sections; 

c) by means of a Kirch h off algorithm for simulation of multiple reflections, usino each of said zero-offset sections to 
generate a new zero-offset section which only contains simulated, multiple reflections* that is, an ensemble of sim- 
ulated multiple sections; 

*s d) by filtering, adjusting amplitude, frequency and phase of each of the so-obtained simulated multiple sections to 
each corresponding zero offset section of step b), so as to generate en ensemble of adjusted multiple section* hav- 
ing amplitude, frequency and phase compatible with the registered multiple reflection; 

e) from each section of the ensemble of the zaro-offeet sections of b), subtracting the corresponding adjusted mul- 
tiple section of d), so as to obtain an ensemble of multiple-free sections. 

so 

[0044] According to a second embodiment as illustrated in FIGURE 5, steps c) and d) are modified in the following 
way: 

c) Btaoking a sorted subset of the ensemble of zero offset sections of step b) so as to generate a single stack sec- 
ss tion (pseudo zero offset section), which is used as input to the Kirchhoff algorithm for simulation of multiple reflec- 
tions, so as to generate a single simulated multiple section; 

d) by filtering, adjusting the amplitude, frequency and phase of the so-generated single simulated multiple section 
to each section of the zero offset sections ensemble of step b) so as to generate an ensemble of adjusted multiple 
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sections. 

[0045] It should be understood that in the first two embodiments of the present invention, the process flow as rep- 
resented by Figures 4 and 5 is applied in such a way that the adjusted multiple which is obtained for a certain common 
s offset section is used only as input for the subtraction step, without any further use in other common offset sections, 
[0046] On the contrary, in the third embodiment of the present method, the process flow is applied in such a way 
that the adjusted multiple which i$ obtained for a certain common offset section is used as input for the nearest common 
offset section. 

[0047] Thus, the third embodiment of the present method involves choosing the common offset section having the 
io shortest offset after the transformation of said common offset section to zero offset. The transformed common offset 
section having the shortest offset is thus considered as the first zero offset section. The reason for the choice of the 
shortest offset Is that It has undergone the least transformation and therefore the highest possible accuracy may be 
obtained by applying the zero offset transformation. 

[0048] According to said third embodiment, as illustrated in FIGURE 6. steps c) and d) are modified in the following 
1$ way: 

c) from the ensemble of zero offset (ZO) sections, sorting the first zero-offset section, whloh ic used as input to the 
Kirchhoff algorithm for simulation of multiple reflections, so as to generate a simulated multiple (SM) section; 

d) by filtering, adjusting th» amplitude, frequency and phaaa of the eo-generated simulated multiple section to the 
20 first zero offset section so as to generate the adjusted multiple (AM) section corresponding to the first common off- 
set section, said adjusted multiple section playing the role of simulated multiple section for the following, second 
common offset section, and successively to all sorted common offset sections. 

[0049] The three embodiments described above consist in the application of the method of the invention to common 
28 offset sections. 

[0050] For those above embodiments, on attaining the final stage of the method, that is, the multiple-free section, 
a further, state-of-the-art step may be applied to the obtained multiple-free, that is, to apply the inverse zero-offset trans- 
formation. Such inverse transformation means the inverse NMO correction, the inverse NMO+ DMO correction or the 
inverse MZO process. This means that the source-receiver configuration is recovered, as it was at the stage of common 
30 offset sections. 

[OOSi] Another, fourth possible embodiment is the application of the present method on stack data. Said embodi- 
ment comprises the following steps, as illustrated in FIGURE 7 attached: 

a) applying a proper zero offset transformation and stacking the obtained data in order to generate a stack section; 
33 b) using the so-generated stack section as an input to the KlrchhofT method for multiple simulation, obtaining a zero 
offset simulated multiple section; 

c) by filtering, adjusting amplitude, frequency and phase of the so-obtained zero offset simulated multiple section, 
so as to generate a multiple section of amplitude, frequency and phase compatible with the registered multiple 
reflections, that is, an adjusted muFtiple section; 
40 d) from the zero-offset section (stack section) of a), subtracting the corresponding adjusted multiple section of c). 
so as to obtain a multiple-free section, 

[0052] For all described embodiments, the transformation to Zero offset can be performed by applying 

(i) NMO correction; 

(ii) NMO + DMO correction; 

(iii) or MZO process 

wherein NMO means Normal Moyeout, DMO means Ojp Moveout and MZO ia Migration to Zero Offset. 

so (0053] The filtering prooeee is carried out according to the usual wsy, that is by generating a matching filter using 
the simulated multiple reflection as input and the observed multiple as desired output. The simulated multiple is found 
at the simulated multiple section while the observed multiple is contained at the zero-offset section. The aim of the fil- 
tering is the matching of the simulated multiple reflections to the observed multiple reflections. 
[0054] The filtering step may be performed In several alternative ways, 3uch as; - 

55 

i) by using a single matching filter for the whole section 

ii) by using a set of matching filters, each filter being applied to a portion of the adjusted multiple section; 

ill) by using a singte matching filter for the whole section, effecting filtering and subtraction at the same step; 
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rv) by using a set of matching filters, each filter being applied to a portion of the adjusted multiple section, effecting 
filtering and subtraction at the same step. 

[0055] The Inventive metnod or muiupis simulation using a Kirchhoff algorithm can be performed either: 

5 

(i) using equations (1 ) to (4). It comprises then a summation of figures collected along a stacking line on input data 
such as the pseudo zero offset section; 

(li) or using equations (5) to (8). Jt comprises then a) migrating input data to depth domain so as to obtain a migrated 
section and b) demigrating the so-obtained migrated section so as to obtain a simulated multiple section, 

to 

[0056] The examples of the present invention have been obtained by applying the method of the Invention in one 
single step. However, the same or different examples may be obtained by applying the two-step approach. 
[0057] The present invention will now be illustrated by the following Examples, which should not be construed as 
limiting the Invention. 

15 

EXAMPLE 1 

[0058] The Example corresponds to the application of the inventive method to a PETROBRAS offshore 2-D seismic 
line recorded at Brazilian continental margin. In the region under investigation, average sea-floor depths are i,soo m, 
20 dips varying up to 5.2 degrees. The main parameters relative to the acquisition geometry and processing are displayed 
in Table 1 below. 



TABLE 1 



3S 



30 



Interval between shots 


26.57 m 


Interval between receivers 


26.67 m 


Interval between CMPs 


26.67 m 


Sampling interval 


2 ms 


Processed record time 


6.0 $ 


Minimum offset 


216 m 


Maximum offset 


4191.83 



[0059] Figure 8 shows a portion of a stacked section of the whole seismic line to which the method has been 
applied. 

[0060] The available data set consisted of an ensemble of common-offset sections, each of which having been sub- 
40 jected to a whole pre-stack processing sequence, which included conventional multiple attenuation using F-K methods. 
Figure 8 shows a stack of 22 common-offset sections with offsets up to 778 m. Examination of the common offset sec- 
tion shows that, in the considered smalJ-offset range, the F-K method has not been able to correctly attenuate, free-sur- 
face multiples. Note, for example, that sea-bottom multiples can be clearly observed at times greater than 3.2 $ and 
strongly interfere with primary reflections. As an example, Figure 9 readily shows a primary reflection being cut by the 
45 sea bottom multiple. This could be an important target event for which a clear image would be needed. Therefore, the . 
successful withdrawal of said residual multiples from the ensemble of small-offset sections would demonstrate the 
excellence of the proposed algorithm as an adequate tool in the suppression of multiples. 

[0061] To carry out the actual multiple suppression from the given data* the following strategy has been employed: 

so a) the input date set was sorted into an ensemble of common-offset section a; 

b) all sections of said ensemble were NMO corrected using the constant velocity of water. The NMO-corrected sec- 
tion obtained from the common-offset section of shortest offset has been chosen to be the best approximation of a 
zero-offset section. Note that this choice is not unique. An alternative for said best zero-offset approximation could 
be, for instance, a stacked section of a number of previously NMO-corrected sections; 

ss c) the Kirchhoff-type summation algorithm is then applied to the above approximated zero-offset section to simulate 
a section of zero-offset multiples, honoring, basically, traveltimes; 

d) by means of an adaptive filter, automatically matching amplitudes, frequency and phase of the above-simulated 
multiples to the corresponding ones of the observed multiples so as to yield a new updated section of simulated 
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multiples* now honoring traveltime and amplitude requirements; 

e) subtracting the previous section of simulated multiples from the input section to produce a new, free-multiple sec- 
tion; 

f) the simulated multiples In section (d) were considered as an initiat simulation of the multiples with respect to the 
s next NMO-corrected sections obtained from the original data. Analogous application of the adaptive filter was car- 
ried out on these multiples to match their amplitudes and phases to the corresponding ones pertaining to the new 
input section. The obtained section of simulated multiples was subtracted from the input section to produce its cor- 
responding free-multiple section; 

9) the process was repeated to all common-offset sections under consideration; 
io h) finally, each common-offset section was subjected to an inverse NMO correction and, hence, prepared for the 
next processing steps such as stacking or pre-stack migration. 

[0062] Note in Figure 1 0 the complete suppression of the sea-bottom multiple. The success of the application of the 
algorithm according to the method of the invention in suppressing the sea-bottom multiple from the target reflection and 
is leading to a very clean image may thus be demonstrated. 

EXAMPLE 2 

[0063] Example 2 corresponds to the application of the method of the invention on stack section extracted from a 
20 3D seismic volume. The main acquisition and processing parameters of the selected line are presented in Table 2 
below. 



TABLE 2 



30 



Interval between shots 


25-00 m 


Interval between receivers 


12.50 m 


Interval between CMPs 


6.25 m 


Sampling interval 


2 ms 


Processed record time 


4.0 9 


Minimum offset 


148.00 m 


Maximum offset 


235,50 m 



[0064] The strategy used in this case is described in the following steps: 

a) stacking the selected data after the application of the NMO correction using the water velocity; 
40 b) using the stack section to simulate multiples by means of the Kirchhoff summation, generating the simulated 
multiple section; 

c) adjusting and subtracting the simulated multiple section from the stack section, so as to generate a new multiple- 
free stack section. 

45 [0065] Figure 11 presents the results generated after each one of the three steps described above. Figures 12A 
and 128 show the zoom images of Figure 11. Figure 12A shows the stack section before the application of the method 
while Figure 12B shows the same stack section after the successful suppression of the multiples resulting from the 
application of the inventive method, 

[0066] The Kirchhoff simulation used in the present invention presents the following main advantages; i) It is target 
SO oriented, that is, it can be applied to a definite constrained region, using rays that illuminate this region only; Ji) it leads 
to very fast and robust algorithms, at least on simple models such as constant velocity; Hi) it does not require spatially 
equal sampling and therefore may be used for 2-D as well as conventional 3-D streamer data; iv) there is no dip limita- 
tion for multiple simulation and v) it does not require an estimate of the source signature. 

[0067] Interestingly. US patent 6,720,606, ooiumn 3 tines 7-14 states that the method outlined therein is satisfactory 
ss for use in regions where the line of profile is aligned substantially in the direction of maximum dip and provided both 
primary reflections and multiple reflections are present and are recorded. This constitutes a limitation of the method 
since in the absence of the required condition, the obtained results will be impaired. The present method does not 
present such limitation since the background theory does not contain restrictive premises so that the method may be 
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applied to simulate multiples in a 3-D environment. 

[0068] Further, contrary to the teaohinge of the etate-of-the-ert technique, as expressed for example in U3 patent 
5,729,506 (see for example column 2 lines 63 -64) and US patent 5.587,965, which state that such technique involves 
an initial estimate of the source wavelet, the present method dispenses with such estimate. 
5 [0069J In the specification and claims of the presen t invention the Applicant has used the term "Kirchhoff algorithm" 
or "Kirchhoff method" by analogy to well known Kirchhoff migration methods. However, the analogy of the simulation 
method developed by the Applicant and the state-of-the-art methods relates only to the summation aspect while all 
other aspects have been specifically developed for the multiple simulation. 

[0070] The present invention treats the problem of multiple simulation as an imaging problem, based on the general 

to principles stated in the article by Hubral et ai, mentioned hereinbefore. 

J0071] The inventive method is expected to work well on any seismic section that can be a reasonably approxima- 
tion of a zero-offset section. Very good results were obtained by applying the algorithm to deep-water, short-offset sec- 
tions. As most available multiple suppression methods fail exactly on the near-offset range, the present scheme can be 
an alternative in this important range. 

75 [0072] Further, the present method may be applied to a migrated section, In this case, the simulated multiple from 
pseudo zero-onset section needs to be subjected to the same migration procedure as the one applied to the data. 
[0073] Thus, the specification and attached Figures demonstrate that the present method allows the use of near- 
offset section, close to the source, this possibility being not described nor suggested in the open literature. The images 
generated by these near-offset sections are better resolved than the state-of-the-art images. 

Claims 

1. A method for the suppression of multiple reflections from a marine seismic data set having registered primary and 
multiple reflections, which process comprises the following steps; 

a) sorting said input data set so as to obtain an ensembfe of common offset sections; 

b) transforming each of said common-offset sections into a zero offset section, so as to obtain an ensemble of 
zero-offset sections; 

c) by means of a Kirchhoff algorithm for simulation of multiple reflections, using each of said obtained zerc-off- 
3o sat sections to generate a new zero-offset section which only contains simulated, multiple reflections, that is, 

an ensemble of simulated multiple sections; 

d) by filtering, adjusting amplitude, frequency and phase of each of the so-obtained simulated multiple sections 
to each corresponding zero offset section of step b), so as to generate an ensemble of adjusted multiple sec- 
tions having amplitude, frequency and phase compatible with the registered multiple reflection; 

a? e) from each section of the ensemble of the zero-offset sections of b). subtracting the corresponding adjusted 

multiple section of d). obtaining an ensemble of multiple-free sections. 

2. A method for the suppression of multiple reflections from a marine seismic data set having registered primary and 
multiple reflections, which process comprises the following steps: 

40 

ajsortlng said input data set so as to obtain an ensemble of common offset sections; 
bjtransforming each of said common-offset sections into a zero offset section, obtaining an ensemble of zero- 
offset sections; 

c) stacking a sorted subset of the ensemble of zero offset sections of step b) so as to generate a single stack 
45 section that is, a single pseudo zero offset section; 

d) using said single pseudo zero offset section as input to the Kirchhoff algorithm for simulation of multiple 
reflections, so as to generate a single simulated multiple section; 

e) by filtering, adjusting the amplitude, frequency and phase of the so-generated single simulated multiple sec- 
tion of d) to each section of the zero offset sections ensemble of step b) so as to generate an ensemble of 

so adjusted multiple sections; 

f) from each section of the ensemble of the zero-offset sections of b). subtracting the corresponding adjusted 
multiple section of e), so as to obtain an ensemble of multiple-free sections. 

3. A method for the suppression of multiple reflections from a marine seismic data set having registered primary and 
ss multiple reflections, which process comprises the following steps: 

aborting said input data set so as to obtain an ensemble of common offset sections; 

b) transforming each of said common-offset section into a rero offset section, so as ta obtain an ensemble of 
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zero-offset sections: 

c) from said ensemble of zero offset sections, sorting the first zero-offset section; 

d) using said first zero offset section as input to the Kirchhoff algorithm for simulation of multiple reflections, so 
as to generate a simulated multiple section; 

5 e)by filtering, adjusting the amplitude, frequency and phase of the so-generated simulated multiple section to 

the first zero offset section so as to generate the adjusted multiple section corresponding to the first common 
offset section, said adjusted multiple section playing the role of simulated multiple section for the following, sec- 
ond common offset section, and successively to all sorted common offset sections; 

f)from each of the zero offset sections of b}, subtracting the corresponding adjusted multiple section of d) so 
*0 as to obtain the corresponding multiple-free section, 

4. A method for the suppression of multiple reflections from a marine seismic data set having registered primary and 
multiple reflections, which process comprises the following steps; 

15 a)applying a proper zero offset transformation and stacking the obtained data in order to generate a stack sec- 

tion, that is, a single pseudo zero offset section; 

bjusing the so-generated single pseudo zero offset section as an input to the Kirch h off method for multiple sim- 
ulation, obtaining a zero offset simulated multiple section; 

c) by filtering, adjusting amplitude, frequency and phase of the so-obtained zero offset simulated multiple eeo- 
70 tion - 60 as to generate a multiple section of amplitude, frequency and phase compatible with the registered 

multiple reflections, that is, an adjusted multiple section; 

d) from the single pseudo zero-offset section of a), subtracting the corresponding adjusted multiple section of 
c). so as to obtain a multiple-free zero offset seotion. 

?s 5. A method according to claims 1 , 2, 3 or 4. wherein the transformation to zero offset Is performed by applying NMO 
correction. 



6. A method according to claims 1 , 2. 3 or 4. wherein the transformation to zero offset is performed by applying NMO 
correction followed by a DMO correction. 

30 

7. a method according to claims 1 ( 2, 3 or 4. wherein the transformation to zero offset is performed by applying MZO 
(Migration to Zero Offset) processes. 

8. A method according to claims 1 , 2 or 3, wherein further an inverse zero-offset transformation is applied to the 
33 obtained multiple-free section. 

9. A method according to claim 8. wherein the inverse zero-offset transformation is performed by applying an inverse 
NMO correction. 

40 10\ A method according to claim 6\ wherein the inverse zero-offset transformation is performed by applying an inverse 
NMO correction followed by a DMO correction. 

11. A method according to claim 8. wherein the inverse zero-offset transformation is performed by applying inverse 
MZO {Migration to Zero-Offset) processes. 

45 

12* A method according to claims 1. 2 ,3 or 4. wherein the Kirchhoff multiple simulation is performed in one etop. 

13, A method according to claim 12. wherein the Kirchhoff multiple simulation in one step comprises a summation of 
figures collected on input data along a stacking line. 

so 

14, A method according to claims 1 , 2 .3 or 4, wherein the Kirchhoff multiple simulation is performed in two steps. 

15. A method according to claim 14. wherein the Kirchhoff multiple simulation performed in two steps comprises; i) 
migrating input data to depth domain so a? to obtain a migrated section an* II) demlgratlng the So-obtained 

55 migrated section using the multiple demigrating approach. 

16. A method according to claims 1, 2. 3 or 4. wherein filtering is performed by using a single matching filter for the 
whole section. 
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17. A method according to claims 1 , 2, 3 or 4, wherein filtering is performed by using a set of matching filters, each filter 
bemg applied to a portion of the adjusted multiple section. 

18. A method according to claims 1, 2. 3 or 4. wherein filtering and subtraction are performed at the same step using 
5 a single matching filter for the whole section. 

10. A method according to claims 1, 2. 3 or 4, wherein filtering and subtraction are performed at the same step using 
a set of matching filters* each filter being applied to a portion of the adjusted multiple section, 

10 
15 
20 
23 

30 

as 

40 
45 
SO 
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FIGURE 5 



C input data C 



Common Offset Sort 



Common Offset Sections 



Transformation to Zero Offset 



Zero Offset Sections 





Stack 









Stack Section 



Kirchhoff Multiple Simulation 



Simulated Multiple Section 



Adaptive Filtering 



Adjusted Multiple Sections 



TTT 



Subtraction 



Multiple-free Sections 



□ 



PAGE 34/42 1 RCVD AT 2/25/2004 2:07:45 PM [Eastern Standard Time] * SVR:USPT0-EFXRF-1/1 * DNIS:8729306 * CSID:+281 589 1482 * DURATION (mm-ss):15-26 



FEB 25 '04 14:21 FROM:PGS CORP 



+281-589-1482 



T-030 P. 35/42 F-703 



EP 1094 339 A2 



FIGURE 6 
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FIGURE 7 



Input data <^ 







Transformation to Zero Offset 






Stack 







StecK Section 







Kirchhoff Multiple Simulation 







c 



Simulated Multiple 
Section 



Adaptive Filtering 



c 



L 



Adjusted Multipl 
Section 



Subtraction 



Multiple-free 
Section 



PAGE 36/42 ' RCVD AT 2^5/2004 2:07:45 PM [Eastern Standard Time] * SVR:USPT0-EFXRF-1/1 ' DNIS:8729306 ' CSID:+281 589 1482 * DURATION (mm-ss):15-26 



FEB 25 '04 14:21 



FROM:PGS CORP 



+281-589-1482 



T-030 P. 37/42 F-703 



BP 1 094 339 A2 




PAGE 37/42 * RCVD AT 2/25/2004 2:07:45 PM [Eastern Standard Time] * SVR:USPTO-EFXRF-1/1 1 DNIS:8729306 * CSID:+281 589 1482 * DURATION (mm-ss): 1 5-26 



FEB 25 '04 14:21 FROM:PCS CORP 



BP 1 094 339 A2 



+281-589-1482 



T-030 P. 38/42 F-703 



FIGURE 9 




PAGE 38(42 ■ RCWD AT 2C5JW04 2:C7:4S PM [Eastern Standard Time] ' SVR:USPT0-EFm-1/1 " DKIS:8ra306- CSIB:t281 S89 1482 ' DURAH0M (mfiV5S):15-26 



FEB 25 '04 14:21 FROMrPGS CORP 



BP 1 094 339 A2 



+281-589-1482 



T-030 P. 39/42 F-703 



• CDP 



FIGURE 10 



I 

*J3 



PAGE 39/42 ' RCVD AT 2/2512004 2:07:45 PM [Eastern Standard Time] * SVR:USPT0-EFXRM/1 » DNIS:8729306 * CS1D:+281 589 1482 » DURATION (mm-ss):15-26 



FEB 25 '04 14:22 



FROM:PCS CORP 



+281-589-1482 



T-030 P. 40/42 F-703 



EP 1 094 339 A2 




PAGE 40/42 k RCVDAT 2/25/20O4 2:07:45 HI [Eastern Standard Time]* SVR:USPT0£FXRM/1 * DNIS:8729306 ■ CSID;+281 589 1482 "DURATION (mm-ss):1M6 



FEB 25 '04 14:22 FROM:PCS CORP 

to AH*' » 



«Bl-589-1482 



T-030 P. 41/42 F-703 



BP 1 094 339 A2 



I time i 

L J 



S 



g 

. i ■ 



. t . 



9 



§ 



V 



§ 



O 



Si 



> 



1 • J 1 1 

g 


g 9 1 


8 


§ 


§ 


k 


i 


§ 


S 3 i 


5 


§ 


i 




s 



[ 

J ° i 




3 



PAGE41/42 ' RCVDAT 2/2512004 2:07:45 PM [Eastern Standard Time] * SVR:USPT0-EFXRF-1/1 * DNIS:8729306 ' CSID:+281 589 1482 * DURATION (mm-ss):15-26 



FEB 25 '04 14:22 FROM:PGS CORP +231-589-1482 T-030 P. 42/42 F-703 



r 




Title: PSEUDO-OFFSET MIGRATION 


Serial No,: 09/876*9 10 


Attorney: E. Eugene Tblgpen 


Filing Date: June 8, 2001 


Docket No. P30290US 


First Named Inventor: WeizhoneWane 


Date: May 6, 2003 



Contents of Package: 



- Transmittal Form 
Transmittal Letter 

- Information Disclosure Statement with cited Reference 

- Return Postcard 

PTO COPY 



Title: PSEUDO-OFFSET MIGRATION 


Serial No.: 09/876,910 


Attorney: E. Eugene Thtenen 


Filing Date: June 8, 2001 


Docket No. P30290US 


First Named Inventor: Weizhong Wang 


Date: May 6, 2003 



Contents of Package; 



.Transmittal Form 
Transmittal Letter 

Information Disclosure Statement with cited Referen 
Return Postcard 




PTO COPY 




